Adiabatic Elimination in Compound Quantum Systems with Feedback 
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Feedback in compound quantum systems is effected by using the output from one sub-system 
("the system") to control the evolution of a second sub-system ("the ancilla") which is reversibly 
coupled to the system. In the limit where the ancilla responds to fluctuations on a much shorter 
time scale than does the system, we show that it can be adiabatically eliminated, yielding a master 
equation for the system alone. This is very significant as it decreases the necessary basis size for 
numerical simulation and allows the effect of the ancilla to be understood more easily. We con- 
sider two types of ancilla: a two-level ancilla (e.g. a two-level atom) and an infinite-level ancilla 
(e.g. an optical mode). For each, we consider two forms of feedback: coherent (for which a quan- 
tum mechanical description of the feedback loop is required) and incoherent (for which a classical 
description is sufficient). We test the master equations we obtain using numerical simulation of 
the full dynamics of the compound system. For the system (a parametric oscillator) and feedback 
(intensity-dependent detuning) we choose, good agreement is found in the limit of heavy damping of 
the ancilla. We discuss the relation of our work to previous work on feedback in compound quantum 
systems, and also to previous work on adiabatic elimination in general. 
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I. INTRODUCTION 

The quantum theory of continuous Markovian feed- 
back is now well understood Continuous feedback 
arises in a situation where a system continuously inter- 
acts with its environment, and the environment is delib- 
erately engineered such that the influence of the system 
on the environment acts back on the system at a later 
time. This can be described as a Markovian process when 
(a) the natural coupling of the system to the environment 
is approximately Markovian, and (b) the effective time- 
delay in the feedback process is negligible compared to 
any relevant time scale of the system. If the Markovian 
approximation is appropriate, this leads to the great sim- 
plification that the system evolution may be described by 
a master equation of the Lindblad form Q . 

It is possible to divide quantum feedback into two cat- 
egories, which we may call coherent and incoherent, fol- 
lowing Lloyd |H (but without being limited by his defi- 
nitions). In the latter case of incoherent feedback, it is 
not necessary to use a quantum description of the en- 
tire feedback loop. Rather, at some point, it is permissi- 
ble to change from a quantum to a classical description 
by invoking a measurement step. In a quantum optical 
context, this corresponds to electro-optical feedback 
where a photocurrent derived from detecting the light 
radiated by the system is used to control electro-optical 
devices which change the behavior of the system. In the 
former case of coherent feedback, a quantum description 
of the entire feedback loop is necessary. In a quantum 
optical context this corresponds to all-optical feedback 
Q in which the light radiated by the system is reflected 
so that in interacts with the system again, perhaps via 
some other system. 

Continuous quantum feedback may be non-Markovian 



for a number of reasons. The coupling to the environment 
may be non-Markovian. The time delay in the feedback 
loop may be non-negligible. The feedback may act via 
a second system, the ancilla. In this paper we are con- 
cerned with the last possibility. This is of interest be- 
cause it arises very naturally in quantum optics in both 
all-optical [Q and electro-optical |7| contexts. In princi- 
ple, this sort of feedback can be described as a Marko- 
vian process in the larger state space of the system plus 
ancilla. In practice, this procedure is often not useful, 
because of the critical word larger in the previous sen- 
tence. If the required basis size of the system and ancilla 
are N and M respectively, then the Liouvillian for the 
compound system has of order N 4 M 4 elements. Clearly 
for M large, this is much larger than a Liouvillian for the 
system alone. 

Consequently it would be an advantage to obtain a 
master equation for the system alone, without the an- 
cilla. This is possible if the ancilla can be adiabatically 
eliminated. That is, if the ancilla has a decay rate much 
faster than any relevant system rate, so that it is always 
in a steady state determined by the system state. It is 
the purpose of this paper to determine numerically the 
conditions under which this is possible, and to derive the 
resultant master equations under those conditions, for a 
variety of general feedback systems. 

Previous work in this area has left the situation some- 
what confused. Wiseman and Milburn Q considered all- 
optical feedback via an ancilla system, and adiabatically 
eliminated the ancilla. This was shown to be equiva- 
lent to electro-optical feedback for quadrature feedback. 
However, for intensity feedback it was the same only to 
second order in the feedback strength. Moreover, the 
master equation derived (to second order) was not of the 
Lindblad form. 
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Slosser and Milburn |7]] considered electro-optic feed- 
back of the photocurrent from the idler mode of a non- 
degenerate parametric oscillator onto the pump mode. 
Here the signal and idler mode formed the system and 
the pump mode was the ancilla. The procedure they 
adopted for deriving a master equation for the system was 
as follows. They expanded the feedback master equation 
for the compound system to first order in the feedback 
strength, adiabatically eliminated the pump mode, but 
the final result presented for the system master equation 
contained first and second order terms. As in Ref. [Q, 
this second-order master equation was not of the Lind- 
blad form. Furthermore, the steady-state field averages 
were calculated using an unstated aZZ-order master equa- 
tion (which was of the Lindblad form). There are other 
problems with this paper but they are not relevant 
to the present work. 

In this work we show how adiabatic elimination can 
be done rigorously in compound quantum feedback sys- 
tems such as those of Refs. As well as being of in- 
terest in the field of quantum feedback, the methods we 
use for adiabatic elimination are of more general interest. 
While adiabatic elimination of an ancilla mode which is 
linearly coupled to the system is well understood, adia- 
batic elimination with a nonlinear (e.g. proportional to 
the intensity) coupling is not. In particular, the meth- 
ods we use here put the results obtained by Dohcrty and 
co-workers [|lOj on the motion of an atom coupled to a 
damped optical cavity mode on a more rigorous footing. 

This paper is organized as follows. In Sec. II we con- 
sider simple direct-detection feedback, and the four types 
of analogous feedback in compound systems: electro- 
optic feedback via a two-level atom, electro-optic feed- 
back via an optical mode, all-optical feedback via a two- 
level atom, and all-optical feedback via a mode. We show 
that in all four cases it is possible to eliminate the ancilla 
under suitable conditions, giving a master equation for 
the system alone. In Sec. Ill we compare the stationary 
state of these master equations with the solution of the 
full dynamics of the compound systems. For this test 
we choose the free dynamics of the system to be that of 
a below-threshold parametric oscillator, the quantity be- 
ing fed back to be the intensity, and the quantity being 
controlled by the feedback to be the detuning. We also 
compare the results of all five feedback mechanisms with 
that caused by an analogous "reversible feedback" gener- 
ated by a x^ 3 ) nonlinearity. In Sec. IV we conclude with 
a discussion of our results, and present a generalization 
of the all-optical case to multiple optical modes. 



II. ADIABATIC ELIMINATION 

A. Simple Feedback 

In order to discern how the dynamics of a system are 
affected by a feedback loop that includes an ancilla, it is 



useful to know the master equation for simple feedback. 
By simple feedback it is meant that the measurement re- 
sults, based on continuous observation of a source system, 
are immediately used to alter the evolution of the source 
without the involvement of any other quantum system. 
To use an example from quantum optics, a photodetec- 
tor may register photon arrivals from a cavity at discrete 
times and, at these times, some specified change to the 
system may be made (see Fig. [I]). Types of changes in- 
clude altering the optical path length or damping rate 
of the cavity. In the remainder of this paper we will 
often use quantum optics terminology, but it should be 
remembered that the theory is not restricted to optical 
physics. 

The most general form of the simple feedback master 
equation has been derived by Wiseman ||. Consider a 
system with Hamiltonian H and some dissipation at rate 
7 and with lowering operator c. With h set equal to 
unity, the master equation is 

p(t) = -i[H,p}+7V[c]p, (2.1) 

where the Lindblad Q superoperator is 

V[c]=J[c]-A[c], (2.2) 

where for arbitrary operators A and B, 

J[A]B = ABA^ ; A[A]B = %{A% B}. (2.3) 

It is the dissipation which allows for continuous observa- 
tion, the result of which is a current I(t). In this paper 
we are concerned with what is known as direct detection 
where 

I(t) = dN(t)/dt, (2.4) 

where dN(t) is the point process (the increment in the 
number of photons counted) defined by 

[dN(t)} 2 = dN(t) 

E[dN(t)] = -ydtTr[Jcp c (t)] (2.5) 

Here E denotes a classically probabilistic expectation 
value, while the c subscript denotes that the state p c 
is conditioned on the previous measurement results. We 
have assumed that the detection is perfectly efficient; the 
generalization to inefficient detectors is trivial Q . 

Simple feedback arises from adding a Hamiltonian to 
the system evolution of the form 

flfb(t) = I{t)Z (2.6) 

where Z is an Hermitian system operator. Taking into 
account the singularity of I(t), and the fact that the feed- 
back must act after the measurement, it is possible to 
derive a master equation for the system with feedback, 
averaging over all realizations of the stochastic measure- 
ment record I{t). The result is 

p= -t{H,p} +1 V[e~* z c}p. (2.7) 
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To compare this master equation with those obtained 
later it is useful to expand the exponentials to third or- 
der 



C[H] P + 1 V[c]p 

+ 1 {c[Z] + \{C[Z]f + \{C[Z]f) J r\j, 



(2.8) 



where CL4]B = — i[A,B] for arbitrary operators A and 
B. 

The derivation outlined above for the feedback mas- 
ter equation treats the photocurrent I(t) as a classical 
stochastic process, which causes the conditioned system 
state p c to undergo stochastic evolution (known as a 
quantum trajectory There is an alternative deriva- 

tion which treats the photocurrent I(t) as an operator. 
This derivation works in the Heisenberg picture, where 
the system evolution is described by stochastic operator 
differential equations known as quantum Langcvin equa- 
tions |Q. This method is useful for adiabatic elimina- 
tion, so we will briefly review its features. 

Quantum Langevin equations (QLE) are constructed 
without using the con cept of measurement. The dissipa- 
tive evolution of Eq. (2.1) can be derived in a quantum 
optical context from a linear coupling (in a rotating frame 
and with the rotating wave approximation) 

V = i^[v\t)c-c i v(t)] (2.9) 

between the system and a bath of harmonic oscillators. 
Here v(t) is the bath annihilation operator at the point 
at which it interacts with the system. Just before this 
point, the bath is an input vacuum, with field operator 
v in (t) satisfying Q 



[v in (t),vl(t')] = S(t-t% 



(2.10) 



and all normally-ordered moments vanishing. Just after 
this point, the bath is an output (non- vacuum) with field 
operator fll2|] 

«aut(*)=Win(*) + V7c(*). (2-11) 

The photocurrent operator I(t) is simply the intensity of 
the output field 



I(t)=vl at (t)v out (t). 



(2.12) 



Adding together the evolution due to H, V, and -fffb, 
and again noting that the feedback must act after the in- 
teraction, one can derive the following quantum Langevin 
equation for an arbitrary system operator s |3j 



ds = h'l + V7c f ](e iZ se- iZ - s)[v in + ^jc]dt 
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+7(c^sc— —sc'c— —c^cs)dt 

-Vl\ dv L c - ctdy in, s] + i[H, s]dt, (2.13) 

where dV m = v m dt. All operators have time argument t. 
When the expectation value of this equation is taken an 
equation is obtained that can be converted to the master 
equation (^jj) for simple feedback. If Z is set to zero 
then the Langevin equation describes damping alone. 



B. Electro-optic Feedback via an Atom 



The simplest possible ancilla system is a two-level atom 
(TLA). In this section we consider incoherent (electro- 
optic) feedback via this ancilla. The output from the 
system is monitored by direct detection, the results of 
which are used to affect the evolution of the two level 
atom which is coupled to the system, as shown in Fig. |^. 
The system and ancilla are assumed to have approxi- 
mately the same resonant frequency. If the atom is to be 
adiabatically eliminated, it must be heavily damped, in 
which case it will mostly be in the ground state. Then 
the most natural form of feedback involves flipping the 
state of the TLA whenever the photodetector monitoring 
the system makes a detection. This can be achieved with 
a feedback Hamiltonian of the form 



Hfb = -<r x I(t). 



(2.14) 



Here a x is the usual Pauli spin matrix for describing an 
atomic state jl3| . It could be realized experimentally by 
very briefly driving the atom with a pulse of on-resonance 
radiation (a V pulse) which will flip it from the ground 
to the excited state. 

With this form of feedback, the obvious coupling of 
the atom to the system to consider is one proportional 
to the excited state population operator er^er. Here 
a = (°x — *°y)/2 is the atomic lowering operator. Specif- 
ically, 



H, 



coupling 



ScrK, 



(2.15) 



where K is an arbitrary Hermitian system operator. 
When feedback onto the atom in the ground state oc- 
curs the upper state population jumps to a value of 1 
and then decays away, due to coupling to the continuum 
of electromagnetic field modes. In other words, a^a will 
tend to follow the photocurrent. Thus there is a strong 
similarity to simple feedback, if K is chosen to be some 
scalar multiple of Z . 

It is not hard to generalize Eq. (2/7) to include the 
TLA ancilla 



W 



-i[H, 



system 



crW, W] 



+V[exp(-i-(T x )c}W + TV[a]W, 



(2.16) 



where L is the damping rate of the atom and W is the 
density matrix for the compound system. The damping 
rate 7 of the system has been set equal to unity without 
loss of generality. Of course, the operators are now act- 
ing in the joint Hilbert space of the two systems so that 



1. 



and a = 1 



system 1 



(J, etc. 
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The above master equation gives the evolution of the 
density operator for the compound system. At any time 
a partial trace of this operator over the atom could be 
performed to obtain the reduced density matrix for the 
system alone. However, in general, this cannot be done 
to the master equation itself in order to obtain a master 
equation for p sys tem(i)- The obvious exception to this is 
the case where K = and the system is unaffected by 
the atom. 

It is logical that a master equation for the system can- 
not be derived if the atom observables fluctuate, in re- 
sponse to the feedback, on the same time scale as the 
system observables. The effect of feedback would then 
depend on the constantly fluctuating state of the atom 
which, in turn, depends on previous feedback. Removing 
the atom operators from the master equation without re- 
moving information concerning the system is impossible 
due to the coupling that exists between them. Of course, 
a non-Markovian expression could be written down for 
the atom in terms of the system, but this would not lead 
to a Lindblad master equation without some further ap- 
proximation. 

If the atom reacts very quickly to the feedback and re- 
turns to its initial state before more feedback arrives (the 
next photodetection) then this well defined behavior can 
be built into a master equation for the system alone. In 
essence, the atom's state is approximated by its equilib- 
rium value with respect to the instantaneous state of the 
system and operators are replaced by their steady state 
expressions. This procedure is known as adiabatic elimi- 
nation of the atom. 

To proceed with the adiabatic elimination it is noted 
that the total density matrix can be expanded as 



-i[H s ,p s ]+V[c]p s -i[K,p 2 }. 



(2.22) 



w= A><8>|i>UI + pi®UXu 
+ P t i®IT>ai + p 2 ®|TXTI, 



(2.17) 



where the ps exist in the system subspace. All possible 
states of the atom have been included (||) and ( M corre- 
spond to the excited and ground state respectively) . This 
approach is particularly appropriate because of the small 
basis involved. If the above expression for W is substi- 
tuted into the master equation then the atom operators 
can act on the states of the atom. If the coefficients of the 
various orthogonal states are equated the following equa- 
tions for the ps are obtained (the subscript 's' indicates 
the system): 



p = C[H s ]p + J[c]p 2 - A[c]p a + Tp 2 , (2.18) 
p\ = C[H s ] Pl + i Pl K + J[c]p\ - A[c] Pl - -p x , (2.19) 
P2 = C[H s ]p 2 + J[c] Po - A[c]p 2 - T P2 . (2.20) 



By tracing Eq. (|2.17| ) over the atom the reduced density 
operator for the system is 



Ps = PO + P2 

and its evolution equation is found to be 



(2.21) 



Without some approximation this is as far as the elim- 
ination of the atom can be taken. It is not a master equa- 
tion due to the dependence upon p 2 . As discussed pre- 
viously, the limit in which the atom returns very quickly 
to the ground state after feedback needs to be consid- 
ered. Because the probability for photodetection in any 
infinitesimal time period scales as the size of the period, 
the atom is in the ground state for almost all time. The 
approximation that p s sa P q is therefore made. To obtain 
a master equation, an ex pression for p 2 in terms of po is 
needed. From Eq. ( 2.20 ) it can be seen that if T is large 
compared to the other co-efficients of p 2 (except possi- 
bly K) then fluctuations in this operator will be quickly 
damped out and p 2 can then be set to zero. The effect 
of K is to cause rotation of p 2 but not to affect its size. 
The physical picture already described is consistent with 
r being large. Assuming K ~ T 3> 1 (where K ~ V 
means that the operator K scales like T), we find the 
steady state of p 2 to be 



p 2 = (T-C[K})- 1 J[c] Po . 



(2.23) 



When this is substituted into Eq. ( [2.22 ) the master equa- 
tion for the system alone is obtained. With Z = K/T it 



p s = {C[H a ] + V[c] + C[Z](1 - e^])- 1 ^]}^ (2.24) 

It is not immediately clear that this master eq uati on is 
of the Lindblad form j^] . However in appendix A 1 it is 
shown that it can be written as 



p s = -i[H s , p s ] + I dqe-iVie-^c}! 
Jo 



(2.25) 



Some feeling for the nature of the master equation 
can be obtained by an expansion to third order in Z 
(a small feedback approximation) . This gives (subscripts 
dropped) 

p~C[H}p + V[c]p 

+ {C[Z] + (C[Z}) 2 + (C[Z}f}j[c}p. (2.26) 

These terms can be compared to the third order expan- 
sion of Eq. ( |2.8D , with 7 = 1. The difference in second 
and higher order terms means that for large feedback the 
two systems will be significantly different. 



C. Electro-optic Feedback via a mode 

The more challenging task of adiabatically eliminating 
an ancilla that has an infinite number of basis states is 
now considered. Optically, this could correspond to a 
single-mode cavity. The method of expanding the com- 
pound density matrix in terms of the lower number states 
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of the ancilla is not appropriate due to the type of feed- 
back that is utilized. Instead we use Quantum Langevin 
equations, which place no such restriction on the excita- 
tion of the ancilla. 

The output field from the system is once again con- 
tinuously monitored using direct detection (see Fig. ||). 
We take the feedback to be linear driving of the ancilla 
cavity. This causes a jump in amplitude of the ancilla 
cavity when there is a photodetection. It is described by 
the feedback Hamiltonian 



H fb = -(-ib- 



(2.27) 



where b is the annihilation operator for the cavity, e rep- 
resents the amplitude of the coherent driving field and 
I(t) is the operator for the photocurrent output from the 
system. Its effect can be determined from the Heisenberg 
equation of motion for b, 



9fb 



-i[b,Hi 



fb 



(2.28) 



Since I(t) consists of S functions, it is clear that the cav- 
ity field amplitude changes by an amount e/2 whenever a 
photodetection occurs. Note that here the implicit equa- 
tion of motion for b is sufficient to determine its evolution 
because the stochastic term is not dependent upon b [|l4| . 
To provide a feedback circuit that is classically equiva- 
lent to simple feedback in the limit of large damping of 
the cavity, the following choice of coupling is made: 



V 



K 



(2.29) 



where dU m = u m dt. The vacuum field input for the 
driven cavity, u- m , has the same properties as v- m . 

To adiabatically eliminate the cavity, in the limit of 
heavy damping, a QLE will fi rst be determined for a sys- 
tem operator, s. Eq. ( 2.31 ) is greatly simplified, as s 
commutes with all driven cavity operators, to give 



ds = V[J]sdt - [dVl 



i[-{b + tf) + H s ,s]dt. 



(2.32) 



From this it is evident that an expression for b is required 
if a master equation for the system alone is to be derived. 
The QLE for b is 

b = -i— - ^-b - Vfu in + l^utZW- (2-33) 

For large T the fluctuations in b due to system opera- 
tors will be quickly damped out. However, the stochastic 
terms have an infinite bandwidth, so that it is not strictly 
possible to slave an operator that only responds to a finite 
bandwidth, T, to these fluctuations. Although this prob- 
lem can be side-stepped Q it will prove advantageous to 
use the following equilibrium value of b 



b=-—- I dre- TT ' 2 



Vfu in (t - t) 



(2.34) 



The eq uivale nce can be seen if linear damping is included 
in Eq. ( 2.28) ) . The slaved value of b (in the limit of large 
damping 6fb is set equal to zero) is then substituted into 
the coupling, which leaves it in the same form as a simple 
feedback Hamiltonian, given an appropriate choice of K. 
The total master equation is 



W = -i[—{b + tf) + H a ,W] 

+V[e<- b+b " )/2 c}W + TV[b}W, 



(2.30) 



where once again W is the density matrix describing the 
compound system and the damping of the system has 
been set equal to unity. The damping rate of the ancilla 
cavity is given by T. The quantum Langevin equation 
that corresponds to this master equation can be found 
by extending Eq. (2.13). The result for an arbitrary op- 
erator r from either sub-system is 



dr = vt ut 



e£ (6-6t)/2 re -e(6-6t)/2 _ f 



Wout dt 



+V[J]rdt - [dV-lc - cUv in ,r] 
+TV[tf]rdt - VT[dlllb - bUU- m , s 

+i[—(b + tf)+H s ,r}dt, 



(2.31) 



The integral serves to determine the present contribu- 
tion to b from the stochastic terms at time t — r. This 
contribution falls off at rate T/2, the amplitude decay 
rate for the ancilla cavity. The term that is not under 
the integral comes from K which is not stochastic and is 
therefore slowly varying compared to the highly damped 
cavity operators. Thus, b can follow its evolution to a 
very good approximation. 

To simplify matters the Langevin equation for s will 
now be rearranged before substitution so that u- m will an- 
nihilate the vacuum when the expectation value is taken. 
This gives 

ds = V[J]sdt - [dV-lc - cUV in , s] 

+ l -(tf[K,s] + [K,s]b)dt + i[H a ,s]dt. (2.35) 

This is valid as b and commute with system operators. 
We cannot move the stochastic part v- m (t) of v ut{t — t) 
through the system commutator term to annihilate on 
the vacuum. However, it is possible to move the pho- 
tocurrent itself at time t — r as it commutes jl2| . If the 
integrals that will annihilate on the vacuum when the 
trace over the bath is taken are ignored, then we are left 
with 
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if r°° 1 

d = 2 [K ' S] I dTe ^ T ' 21 ^ T ) - 2f [K > [K > S]] 
+£>[c t ]s - [vie - Jv m , s] + i[H s , s}. (2.36) 

If the limit r — > oo is taken the integral reduces to 
2I(t)/T. The resultant equation for s is an implicit equa- 
tion as it was derived by idealizing the properties of the 
cavity and environment fjlf| . An explicit equation is now 
required. 



The term that needs to be treated in Eq. (2.3C) can be 
written as 



^implicit 



eIC[K]s 
f ■ 



(2.37) 



This gives an explicit increment of the form p,[l5| 

rfscxpiicit = dN[exp(-eC[K]/T) - l]a, (2.38) 



where dN = Idt = dN 2 = vl ut v out dt. Remembering that 
the photocurrent is actually evaluated at a slightly earlier 
time than the system operators allows v ou t to be moved 
to the right of the expression. If we put Z = eK/T, 
in order that our equations can be compared to simple 
feedback, then the total Langevin equation is 

ds = [v} n + J\{e iZ se~ iZ - s)[v m + c]dt 



\v-; n + c ue s. 

r 

-[dV^c - cUVin, s] + i[H s , s]dt. 



[Z, [Z,s]]dt + V[J]sdt 



(2.39) 



When the expectation value is taken the stochastic part 
annihilates on the vacuum and the following master equa- 
tion is obtained 

p = -i[H s , p] + V[e- lZ c]p + ^V[Z]p. (2.40) 

The only difference from simple feedback is the third 
term. This is a term of second order in the feedback 
operator Z, and represents a type of noise that will tend 
to smooth over the interesting behavior of the system. 
Clearly it can be made arbitrarily small if e is made large 
enough. A more detailed discussion of this term is given 
in Sec. ElIC 



D. All-optical Feedback via an atom 

We turn now to coherent, or all-optical feedback. Once 
again we begin with the simplest possible ancilla, a two- 
level atom. All-optical feedback via an atom involves 
the reflection of the output field from the system onto 
the atom, where the atom is reversibly coupled to the 
system. Here, the resonant frequencies of the two sys- 
tems are taken to be equal. It is different from electro- 
optic feedback as there is no measurement ste; the light 
is just reflected around a loop with the use of mirrors 
(see Fig. ||) . The theoretical description of such systems 



was developed largely by Carmichael [[L6| and Gardiner 
JItJ and has been termed Cascaded Open Systems the- 
ory. If linear bath-system couplings are assumed then 
the compound master equation is 



W = - i[H s + V, W] + V[c]W + YV[a]W 
+ VT([cW,o-i] + [a, WJ]) . 



(2.41) 



The system damping has been set equal to unity as usual 
and r is the damping rate of the atom. 

In order to investigate the degree to which all-optical 
feedback can replicate electro-optical simple feedback, a 
coupling is chosen that is linear in the excited state pop- 
ulation of the atom. We expect this operator to follow 
the output photocurrent from the system. That is, we 
assume a coupling 



V = Ka^a 



(2.42) 



ident ical to that in Sec. |IIB| . Making the expansion of 
Eq. fl2.17p gives the following for the ps 

po = C[H s ]p + V[c\p + T P2 + VT(c Pl + prct), (2.43) 
p\ = C[H s ] Pl + V[c] Pl + V? (p 2 - p )c f 



+ ipiK- -pi, 

p 2 = C[H s ]p 2 + V[c]p 2 - Vf (epi + ( o :L c t ) 
- i[K,p 2 ) - Tp 2 . 



(2.44) 
(2.45) 



The above equations lead to an equation of motion for 
the system density operator of 



p = C[H s \p + V[c]p-i[K,p 2 ] ) 



(2.46) 



which is the same as Eq. (2.22). To find an expression 
for p 2 the normal procedure of taking T large compared 
to C[H S ] is performed. Thus, pi can be slaved to system 
operators, po an d p 2 - Now as we only require a master 
equation which gives the leading order effect in T _1 of the 
ancilla on the system, p 2 can be set equal to zero in the 
pi equation, which is the approximation po « p. This is 
valid as p 2 ~ po/E By substituting the slaved expression 
for pi into that for p 2 we find after simplification 



P-2 = -yJ 



2iK 



Po- 



(2.47) 



This can now be substituted into Eq. (2.46) to obtain a 



master equation. Writing Z — AK /T, we have 



p = C[H s ]p + V[c]p + C[Z]J 



Zi 



p, (2.48) 



which is the same as the simple feedback Eq. (2 
second order. The third order term is 

\c[Z)(J[Z]-2A[Z])J[c]p. 



to 



(2.49) 
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Again it is not obvious that Eq. (2.48) is i n the Lind- 



blad form, but it is shown in Appendix AS that it can 
be written as 



p = -i[H s ,p]+V 



exp 



Z 

-2i arctan — I <■ 



p. (2.50) 



E. All-optical Feedback via a mode 

The final compound system that will be considered 
involves the output field from a system being reflected 
onto an optical cavity that is coupled back to the system 
(see Fig. |). A Faraday Isolator (comprised of a Fara- 
day Rotator and a Polarization dependent Beam Split- 
ter) prevents reflected light from the cavity returning to 
the system. The only difference in the total master equa- 
tion from the previous section is the replacement of the 
atom lowering operator a with the annihilation operator 
b. Thus a coupling of the form V — Kb% is considered. 

The derivation of a master equation for the system 
alone follows similar lines to that of Sec. II C . The QLE 
for an arbitrary operator is ]f7[ 

dr = +i[H s + V, r]dt + V[c)]rdt - [dV^c - cUV in , r] 
+TV[tf}rdt- Vf[dV^b-bUVi a ,r} 
+Vr (bVc + c)rb - r6 f c - Jbr)dt. (2.51) 

For a system operator this becomes 

ds = V[J]sdt - [dVlc - cUv in , s] 

+ i[H s + Ktfb,s]dt. (2.52) 

The next step is to find an equation for b. The QLE that 
governs it is 



db 



y 2 



Vfv in + Vfc + iKb)dt. 



(2.53) 



This justifies our initial presumption that the cavity pho- 
ton number would follow the photocurrent. For V large 
it is possible to slave b to the system operators and to 
form an integral expressio n for the contribution from the 
stochastic term, as in Sec. II C . The result is 



2 / 2iK 

4 = -7r( 1 + - 

+VT dre- r ^ 1+2lK / r ^ 2 v in (t'T). (2.54) 
Jo 

The same trick of rearranging the QLE for the system 
operator is again used so that, in this case, all of the 
integral terms annihilate. We put 



i[Kb%, s]dt = itf[K, sjbdt. 



(2.55) 



Substituting into this the expression for b and gives 
four terms, only one of which is non-zero when the trace 
over the bath is taken. This term is 



4ic t 



i - 



2iK 



[K,s](l 



2iK 



(2.56) 



In effect, an implicit equation has been derived that has 
no contribution from stochastic operators, resulting in 
there being no need for an implicit/explicit distinction. 
It is now possible to turn the equation for ds into a mas- 
ter equation for the system . Wh en this is done we arrive 
at the same result as Eq. ( [2.48 ). The conclusion is that 
to first order in L" 1 , the cavity has the same effect on 
the system that the atom does, when included in an all- 
optical feedback loop. 

In hindsight, this is what we should have expected, 
as in the limit of large damping only the lowest num- 
ber states of the cavity will be occupied with significant 
probability. One could therefore have expanded the total 
density matrix analogously to the TLA system to obtain 
the same equations immediately. 

The reason why electro-optic feedback onto an atom 
and a cavity were not equivalent is due to the more singu- 
lar nature of the driving of the ancilla. When a detection 
on the system is made the field amplitude of the cavity 
jumps, leading to occupation of higher photon number 
states. These states are, therefore, essential to the de- 
scription of the compound system. Electro-optic feed- 
back onto the atom cannot replicate this behavior. 



III. COMPARISON WITH EXACT RESULTS 



We have shown that in principle it is possible to con- 
sider a variety of different sorts of feedback in compound 
quantum systems, and to adiabatically eliminate the an- 
cillary system to arrive at master equations for the sys- 
tem of interest alone. These master equations should be 
exact in the limit that the ancilla is damped infinitely 
faster than the system. In practice, this will never be 
the case, so it is an interesting question to find out un- 
der what conditions the equations are valid. This can be 
done by simulating the full master equation for the com- 
pound system and comparing to the results of the master 
equation for the system alone. 

To make such a comparison requires specifying the 
feedback operator, Z, and the system Hamiltonian, H s . 
Once this is done, a comparison can be made by look- 
ing at the stationary solutions of the respective master 
equations. While this could be criticised as not being 
a complete test, it has the advantages of definiteness 
and ease of calculation (in some cases at least). Further- 
more, we choose a system (a damped optical mode) and 
Hamiltonians H s and Z such that the stationary solutions 
have enough structure for the comparison to be interest- 
ing. The comparison is both quantitative and qualitative, 
with the use of the Bures distance yl| as a measure of 
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the difference between the state matrices and the Wigner 
function [Jl9| to illustrate them. 

In the hope of getting some interesting states we take 
the system to be a damped single mode optical cavity. 
That is, we choose c = a, an annihilation operator satis- 
fying [a, a'] = 1. We choose a system Hamiltonian (in a 
rotating frame) of 



H 



i\ f 2 



(at) 



h2i 



(3.1) 



This describes a degenerate parametric amplifier ("two 
photon" driving), which can be realized by driving an 
intracavity crystal with a x^ non-linearity with light at 
twice the resonant frequency. For A positive, this results 
in squeezing of the X2 quadrature of the field inside the 
cavity, and stretching of the X\ quadrature. The two 
quadratures are defined in this paper as 



X x = a 



X 2 = —i(a — a') 



(3-2) 
(3.3) 



Without feedback, the master equation with two-photon 
driving and damping will have a stationary solution only 
for A < 1. That is, A is the threshold parameter. 
The feedback operator is chosen to be 



Z 



X a a - 



(3.4) 



We can get a feel for the effect of this type of feedback 
by u sing Z in the simple feedback Hamiltonian given in 
Eq. (Eg). That is, 



Hib(t) = X I{tWa. 



(3.5) 



This represents a detuning of the system cavity propor- 
tional to the photocurrent. It will cause the master equa- 
tion to have a stationary solution regardless of A, as will 
be shown. As the mean photocurrent is equal to the ex- 
pectation value of the photon number operator for the 
system, this Hamilto nian is akin to a x Kerr non- 
linearity p9|. In Sec. Ill F a comparison of feedback to 



such a nonlinearity is made. 



(n\J (e"**) 



iX\ a a 



p\m) 



(3.8) 



Now this particular system has the property that p nm = 
for \n — m\ odd as the two photon driving is the only 
source of coherences. These coherences exist between el- 
ements with I to — m\ even. Hence, if X — qir, with q 
an integer, then the feedback has no effect. Investiga- 
tion into the states produced with a value of feedback 
close to this revealed that they are extremely sensitive 
to any parameter variation. This implies that it is not a 
suitable regime for the testing of adiabatic elimination. 
The most obvious alternative is to choose the maximum 
feedback regime. It is clear that this is achieved with 
X = {q + 1/2) 71 "- The states produced are much less 
sensitive and also have the advantage that, for simple 
feedback, there is no threshold to the driving strength 
above which the photon number becomes infinite. For 
the remainder of the paper we choose x = 7r/2. 

The two-photon driving strength A was chosen to be 
as large as possible, given the constraints on the maxi- 
mum basis size that could be simulated. This amplified 
the interesting effects of feedback. Not surprisingly, the 
simulations of the compound systems are the most com- 
putationally intensive and provide the upper basis size. 
It was found that the limit for the system cavity basis 
size required that photon numbers above 35 had to be 
truncated. For an accurate simulation |pof this gives a 
maximum driving strength of about A = 2.2. Where 
possible, the compound systems were examined in the 
same regime as simple feedback, but for some the driv- 
ing threshold of A = 1 remains in force, so A = 0.97 was 
then chosen. 

The numerical simulations were greatly aided by the 
use of the Quantum Optics toolbox for Matlab As 
noted above, we gauged whether the adiabatic elimina- 
tion is valid by investigating the steady states of the 
systems. The simple feedback system involved a small 
enough Liouvillian that matrix inversion methods can be 
used. The Wigner function of the steady state density 
matrix for simple feedback, with A = 2.2 and x — tt/2, is 
shown in Fig. a. A plot with A = 0.97 is also included. 



A. Simple Feedback 



The master equation for simple feedback is now 



p=~[a 2 -tf)* t p]+V[e 



a\p. 



To simplify the numerical analysis we choose a single 
feedback strength for which simulations will be run. To 
aid this decision the effect of feedback is analyzed. Con- 
sider the following quantity: 



J 



If this is evaluated in the number basis then we get 



B. Electro-optic Feedback via an Atom 

Electro-optic feedback via an atom can be compared 
(3-6) to t he sim ple feedback just considered if we insert in 
Eq. (|2.15[) K = TZ = ga^a, where g = T-k/2. To test 
the adiabatic elimination simulations were run for various 
values of r. It is only for large T that correspondence be- 
tween the full dynamics and the adiabatically eliminated 
master equation is expected. A physical realization of 
this coupling is a far detuned atom in the standing wave 
(3.7) of a single mode cavity |n|. This also introduces a term 
into the system Hamiltonian of the form Sa^a, where S 
is the difference in resonant frequency of the atom and 
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system. It is of interest to determine whether the same 
results are obtained if the adiabatic elimination is done 
at the same time, rather than after, the large-detuning 
approximation is made. This is addressed in appendix [Bj, 
and the answer is affirmative. 

The full master equation is of the form 



W = -i 



iX 



{a 2 - (a r ) 2 } + gcrWa + 5a^a, W 



+V[exp(-i-a x )a]W + TV[a]W. 



(3.9) 



The reduced density matrix for the system at steady state 
needs to be found. Once again, the Liouvillian is small 
enough that we can set W — and solve the equation 
CW = for the non-trivial solution. Simulations were 
run for values of T from 1 to 100, with g altered ac- 
cordingly. Note that the detuning actually has no effect 
on the system dynamics. The reduced density ma trice s 
produced are compared with those found from Eq. (2.24) 
with the aid of the Bures distance, which gives a measure 
of how distinguishable two mixed states (pi and p^) are. 
The Bures distance is defined as [E2| 



3 (pi,p 2 ) = \/2(l-Tr 



P\P2\J P\ 



(3.10) 



All pairs of density matrices of the same size have a Bu- 
res measure that is mapped onto the real numbers be- 
tween zero and Fig. shows how the state produced 
by the compound master equation approaches that pro- 
duced by the adiabatically eliminated master equation. 
As T is increased the Bures distance decreases and the 
Wigner functions become more similar to the adiabatic 
state. This shows that the adiabatic elimination is valid 
in this system for surprisingly small values of T. 

A comparison of the stationary Wigner functions pro- 
duced here with those of simple feedback reveals that 
there exists vast differences between these feedback 
schemes. This is not surprising as it is only to first order 
in Z that the equations are the same, and the param- 
eters we have chosen correspond to Z quite large. The 
most obvious visual differences include the presence of a 
shearing effect and the loss of reflective symmetry in the 
X2 quadrature. 



C. Electro-optic Feedback via a Mode 



In Sec. |I Q electro-optic feedback via a mode was con- 
sidered. In the limit of the ancilla mode being damped 
on a time scale small compared to those of the system, 
Eq. (2.40) was derived. The feedback operator was set 
as Z = eK/T so that we could make a comparison to 
simple feedback. It follows that the system coupling op- 
erator, K, is of the same form as the previous section: 
K = go) a. The coupling V = ga^a(b + b^)/2 could be 
physically achieved via a four wave mixing process in a 



material [jl9|. The fourth field would have to have 
the same frequency as the ancilla cavity for conservation 
of energy. 

No w th at Z has been specified, the third term in 
Eq. ( 2. 40] ) can be discussed more explicitly. This can 
be done by considering the evolution of the phase op- 
erator, which has an approximate commutation relation 
with the number operator of n] = —i p3|j . It can then 
be shown that this term causes phase diffusion at a con- 
stant rate, implying that the features of the state which 
are dependent upon a distinct phase are lost. With the 
notable exception that the photon number is not directly 
affected, there are many similarities with damping. 

For simulation, parameters are chosen so that eg/T = 
tt/2, A = 2.2 and r/2e 2 = 0.001. The last equality main- 
tains the phase diffusion term at a small and constant 
level. This ensures that the same state is always pro- 
duced by the adiabatically eliminated master equation. 

It is worth mentioning how the full dynamics were sim- 
ulated. Due to the jump in the field amplitude of the an- 
cilla cavity when a detection on the output of the system 
is made, the basis size required for an accurate simula- 
tion is large. If the amplitude jumps by an amount e/2 
then the photon number will increase by (presuming the 
initial field was small) e 2 /4. A second detection on the 
system occurring very soon after the first, would push the 
photon number even higher. In fact the computational 
resources available were not su fficient to allow e ven a 

25) of Eq. (fOol). 



quantum trajectory simulation [jll| 
The solution was to make a unitary transformation to a 
frame in which the evolution of the driven cavity due to 
feedback was separated from that due to quantum noise. 
That is, the mean amplitude of the field was described 
classically while the quantum representation of the noise 
was maintained. The unitary transformation used was 



t/ = exp[ e /(f)(6-&t)/ 2 ], (3.11) 
where f(t) is defined by 

/(<)=/ dsexp[-r(t-s)/2]I(s). (3.12) 
J — 00 

Here, I(t) is the c-number stochastic photocurrent. The 
price of a reduced basis size is a time de pende nt Liouvil- 
lian. When the transformation of Eq. (3.11) is applied 
to the implicit master equation (feedback is described by 
a feedback Hamiltonian instead of the exponentials) an 
equation is obtained that is already of an explicit form 
(see appendix |^) 



W 



ga)a{b + + ef{t)} - ^ {a 2 - (a^) 2 } , W 



+T>[a]W + TV[b]W. 



(3.13) 



It can be seen that ef(t) represents the amplitude of 
the driven cavity. Although f(t) is stochastic, it is a 
smoothed (non-singular) version of the photocurrent and 
can therefore be treated without worrying about the 
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stochastic calculus. Note also that since U contains only 
ancilla operators, the system state matrix p — Tr^IF] is 
the same as before, Tr(,[W]. 

The transformed master equation was simulated using 
quantum trajectory methods. It is shown in Fig. |^ that 
as r becomes large the adiabatically eliminated master 
equation becomes a very good approximation to the full 
dynamics. Clearly, though, T has to be pushed to much 
higher levels than the TLA damping for this correspon- 
dence to hold. One reason for this is that the Wigner 
functions of the steady state density matrices for electro- 
optic feedback onto a mode have much greater struc- 
ture, meaning that a measure such as the Bures distance 
(which measures the distinguishability of states) will be 
more sensitive to small differences. It also is likely that 
the parameter regime chosen is one in which this system 
varies quickly, with the result that adiabatic elimination 
will only be valid at very large T. 

A comparison of the Wigner functions [Figs, ^(a) and 
(c)] with that produced with simple feedback [Fig. ||(a)] 
shows the expected similarity. It is expected because 
Eq. ( [2.40 ) only differs from simple feedback due to the 
presence of the double commutator noise term, which was 
chosen to be small. 



E. All-optical Feedback via a Mode 



It was shown in Sec. II E that in the adiabatic limit 



all-optical feedback onto a mode has the same effect as 
feeding back onto a TLA. Therefore, the same threshold 
for the driving strength exists for this system (A = 1). 

The basis size required here is not as large as for 
electro-optic feedback because the photons leak out of 
the system and into the ancilla cavity, giving a smooth 
variation of photon number. Despite this, a quantum tra- 
jectory simulation was still found to be necessary. The 
results obtained for A = 0.97 and ig/T = tt/2 can be 
found in Fig. |l^. The adiabatic state is, of course, the 
same as for all-optical feedback onto a TLA. There is 
a notable difference in the speed at which the full dy- 
namics approaches this state. At low damping the Bures 
distance is already very low. The conclusion is that the 
ancilla mode has minimal effect on the system when in- 
cluded in an all-optical feedback loop. 



F. Comparison with "Reversible Feedback" 
Generated by a x' 3 ' Non-linearity 



D. All-optical Feedback onto an Atom 

The basis size of the TLA ensures that s imul ating the 
full dynamics of all-optical feedback [Eq. (2.41)] is rela- 
tively easy. However, a threshold driving strength exists 
(A = 1) for this system which means that the adiabati- 
cally eliminated master equation cannot be tested in the 
same regime as the previous sections. Instead we set 
A = 0.97 which enabled us to perform an accurate simu- 
lation with the computatio nal re sources available. 

The Hamiltonian of Eq. ( 2.41 ) includes the paramet- 
ric amplifier driving and also the coupling of Eq. ( [2.42 ). 
Once again we choose K = ga^a and set Ag/T = tt/2, 
while varying T and g. T he B ures distan ce be tween the 
states produced by Eq. ( 2.41 ) and Eq. ( 2.48| ) is shown 
in Fig. ^, as are some Wigner functions for the full dy- 
namics and the adiabatic state. It can be seen that the 
state produced with the full dynamics approaches the 
adiabatic state at a similar rate, as T is increased, to 
electro-optic feedback via a TLA. 

There is a large similarity between the state produced 
via simple feedback in Fig. ^ (c) and that in Fig. ^| (a), 
with the presence of shearing being the most notable dif- 
ference. This closer correspondence to simple feedback 
than that of the electro-optic feedback systems is not sur- 
prising given that the adiabatic all-optical master equa- 
tion was the same as simple feedback to a higher order 
(second). The smaller driving also contributes to the 
closeness of the states. 



Finally, we consider the effect of placing a x^ material 
inside an optical cavity driven by a parametric oscillator. 
There is no feedback loop involved. The Hamiltonian 
generated by the x^ non-linearity (a Kerr non-linearity) 
is given by [EM 



H KclI = 



(3.14) 



The Heisenberg equation of motion of the annihilation 
operator due to this Hamiltonian is found to be 



d = —ixia) a)a. 



(3.15) 



Thus, it is clear that the x^ 3 ' non-linearity causes a de- 
tuning proportional to the intensity of the field inside the 
cavity. In this way the system has a self-awareness that 
is similar to simple feedback, which is why a compari- 
son is relevant. In fact, it can be shown that the two 
systems are classically equivalent given the same choice 
of the parameter x- For large feedback the two systems 
diverge when treated quantum mechanically. One of the 
main reasons behind this is that the Kerr effect displays 
no periodic dependence upon its magnitude, whereas the 
simple feedback does. This is illustrated in Fig. |Tl](c) , 
where the Bures distance between the steady states of the 
two systems is plotted for varying \. The Wigner func- 
tion of the "reversible feedback" steady state for x = 7r /2 
with A = 2.2 and A = 0.97 is given in Figs. |ll](a) and (b) 
respectively. They are seen to be very different from any 
of the steady states produced by feedback. 



10 



IV. DISCUSSION 



A. Summary 

The description of feedback in compound quantum sys- 
tems (where the output from the system is used to control 
the evolution of the ancilla, which is reversibly coupled 
to the system) is greatly simplified if the ancilla can be 
adiabatically eliminated. We have shown how this can be 
done for four generic cases, arising from considering two 
forms of feedback (all optical or coherent, and electro- 
optical or incoherent) and two types of ancilla (a two- 
level atom, and an optical mode). The four resulting 
master equations for the system alone are given below. 
They are the most important results of this paper. We 
also include the perturbative expansions of these master 
equations to third order in the feedback operator Z. All 
of the equations are identical to first order in Z, but differ 
in second or third order. 

For comparison, we begin with simple feedback (that 
is, with no ancilla) based on detection of the intensity 
fo!,,* fonnt. of the output field b n nt = b m + c, and using 



out "out of the output field b 

out 



the feedback Hamiltonian 

Ha, = I(t)Z. 
The master equation for this is 



p=-i[H lP \+V[e- lZ 



-i[H,p} + V[c]p + C[Z]J[c}p 



\(C[Z]f 



\(c\zYf\j\c,, 



(4.1) 



(4.2) 



(4.3) 



Here C[Z]B = —i[Z, B] as before. The remaining master 
equations result from trying to reproduce this form of 
feedback via an ancilla. 

The first master equation derived using adiabatic elim- 
ination is for electro-optic feedback via the inversion of a 
two-level atom: 



POO 

-i[H,p]+ / dqe~ q V[e- lqZ c] Pl 
Jo 

-i[H,p} + V[c]p + C[Z]J[c}p 
+ {{C[Z]f + {C[Z]f}j[c] P . 



(4.4) 



(4.5) 



This differs from the simple feedback master equation 
(4.2) at second order in Z. The second is for electro- 
optic feedback via one quadrature of an optical mode: 



f > -- -i[H, p] + V[e- lZ c}p + -^V[Z). (4.6) 

( - 



The expansion of the above equation can be found from 
that of the simple feedback. The size of the extra second- 
order term is determined by T, the damping rate for the 
ancilla mode, and e, the strength of driving of the ancilla 
mode. 



Turning now to all-optical feedback, we have found 
that the same master equation arises regardless of 
whether the feedback is via the inversion of two-level 
atom or the intensity of an optical mode. It is: 



-i[H,p]+V 



exp 



Z 



-2i arctan — | c 
2 



i[H,p]+V[c]p + C[Z]J[c]p 
1 
4 



\{C[Z]) 2 



P, (4-7) 
C[Z]{J[Z)-2A[Z])\j[c)p. 



(4.8) 



Unlike Eq. fl4.4|) , this differs from Eq. ( fl2| ) only at third 
order in Z . 



B. Relation to Previous Work 

As mentioned in the introduction Slosser and Milburn 
perform adiabatic elimination of the pump mode of a 
non-degenerate parametric oscillator. In their system the 
pump mode is driven by the output photocurrent from 
the idler mode. The procedure they adopt is similar to 
that contained in 



Sec. [IB of this paper, in that they 



expand the density matrix in terms of the lower number 
states of the pump mode. However, in Sec. II C we have 



already noted that this is not appropriate when dealing 
with direct detection feedback onto a mode. Higher num- 
ber states are essential to the description of the system 
if the feedback strength is large. For this reason they 
limit the feedback strength to small and moderate val- 
ues, with a generalization to larger feedback contained 
in their appendix. This appendix does not explain the 
origin of the all-orders feedback term. The techniques of 
adiabatic elimination using QLE's that are presented in 
this paper make it easy to treat their system rigorously 
to all-orders in the feedback strength. The final result, 
using their definitions and our superoperators, is (with 
perfect detection assumed): 

p = ela)^ - ab, p] + 2TV[ab]p + "faV[a]p 

+7 b X>[exp ( X ab - X ^tf) b]p. (4.9) 

Note that the second term here is the one analogous to 
the final term in our Eq. ( |4.6| ). 

Doherty and co-workers consider a strongly interact- 
ing system comprised of an atom inside a cavity [fiol . 
The methods used for adiabatic elimination are similar 
to those used in this paper. They form QLE's for oper- 
ators from any of the three subsystems (center of mass 
motion, internal state and the cavity mode) and then set 
the time derivatives of, first, the internal state operator 
and, second, the cavity operator, to zero. They then 
substitute into the QLE for the momentum operator p x . 
After a conversion to the explicit form of the QLE, they 
show that the QLE they derive is compatible with the 
master equation (using their notation) 
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2m 



■P 



-V 



exp I —2i arctan ■ 



P- 



where 



Z = 



5q cos 2 kzx 
Ak 



(4.10) 



(4.11) 



Note the similarity with our equation resulting from adi- 
abatic elimination of an optical mode where the coupling 
is via the intensity, (but of course there is no feedback 
here so our operator c is replaced by the c- number a). 
The derivation of this master equation in Ref. |l(| is not 
completely rigorous in that other master equations would 
also be compatible with the QLE they derive for p x . How- 
ever, it would be straightforward, using the technique we 
introduced in Sec. HE, to make it rigorous. 



The work done on all-optical feedback in this paper 
follows on from that done by Wiseman and Milburn [Q . 
They were able to show that all-optical feedback onto a 
mode could replicate electro-optic homodyne-detection 
feedback, but they could only prove equivalence with 
direct-detection feedback to second order. Here, we have 
shown that this is because the equivalence only holds to 
second order. We have done this by finding the mas- 
ter equation to all-orders in the feedback strength, and 
showing it to be of the Lindblad form. 

Showing that all-optical feedback via an ancilla (be 
it a two-level atom or a mode) cannot replicate electro- 
optical direct detection feedback, leads naturally to the 
question of whether a more complicated all-optical feed- 
back scheme can replicate direct electro-optic feedback. 
Since the feedback is replicated to second order, a fruitful 
approach would seem to be to make the feedback weak, 
while multiplying up the number of ancillae to compen- 
sate. In appendix [d] we consider the case of N ancillae, 
with coupling to the system scaling as 1/N, where the 
output of the system is fed sequentially into all of the 
ancillae. We show that in the limit N — > oo, this hypo- 
thetical all-optical feedback scheme does indeed produce 



the simple electro-optic feedback master equation (4.2 



C. Conclusion 

We have shown that it is possible to greatly simplify 
the description of feedback in compound quantum sys- 
tems by adiabatically eliminating the ancilla, to give mas- 
ter equations for the system alone. In essence, we have 
found the first order in effect of the ancilla upon the 
system, where T is the ancilla decay rate. We have done 
this for a variety of ancillae and forms of feedback, and 
found good agreement with numerical simulations of the 
dynamics for the full compound quantum system. The 
master equations in the various cases are quite different, 
and their range of validity (that is, how large T has to 
be for them to be valid) was also found numerically to 



differ. For the numerical simulations we of course used a 
particular system, but the equations we derive are very 
general. 

The primary motivation for this work is the reduction 
of basis size that is necessary to describe the evolution 
of the system. It is hoped that the derived equations 
will prove to be helpful to co-workers. However, we note 
that numerical testing (to find the regime in which these 
equations are a good approximation) may be necessary 
to determine when it is appropriate to use them. Apart 
from these practical advances, we feel that the previously 
existing confusion in the literature, as discussed in the in- 
troduction, has been resolved, and the procedure of adi- 
abatic elimination in compound quantum systems with 
feedback is now on stable ground. 
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APPENDIX A: PROOF OF LINDBLAD FORM 



1. Electro-optic Feedback onto a TLA 



To show that Eq. ( 2.24 ) can be written in the Lindblad 
form the following identity will first be established: 



(T-C[K]) / dxj 
Jo 



-x(r+2iK)/2 



p = p. 



(Al) 



Multiplying the equation through by two arbitrary eigen- 
states of K, (a\ and from the left and right respec- 
tively, the following is obtained: 



poo 

p a0 [T + i(a - /?)] / dxe- x ^ a ^ = p afi 
Jo 



(A2) 



After the simple integration is performed the identity is 
proved. Before using this the following rearrangement is 
made: 



C[Z]{\-C[Z])-^{\-C[Z])-'-l. 



(A3) 



Upon use of th e identity with Z = K/T the master equa- 
tion Eq. fl2.24j ) becomes Eq. (|2.25|) . 



2. All-optical Feedback onto an Atom 

In order to show that the master equation can be writ- 
ten as in Eq. ( 2.50| ) it is sufficient to show that 



exp [— 2i arctan (Z/2)] p exp [2i arctan (Z/2)] 
= C[Z]J \{X + iZ/2)~ 



(A4) 
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Note that J[c\ has been omitted as it is a multiplica- 
tive factor on both of the superoperators. Consider the 
following non-Hermitian operator that can be put into a 
modulus and argument form: 



1 



1 + (Z/2Y 



(A5) 



That r = 1 can be quickly verified. The argument is 
given by 



9 = arctan (Z/2) . 



(A6) 



By taking the logarithm of |A5| an alternative expression 
for the argument is obtained 



(A7) 



If the logarithmic form of arctan (Z/2) is used, then the 
exponentials of Eq. (A.4) disappear. The LHS of that 
equation becomes 




1 + [Z/2) 2 jl + iZ/2) 2 _ 
(l + iZ/2f P \ + {Z/2f P ' 



(A8) 



It is now noted that 1 + {Z/2) = (1 + iZ/2) (1 - iZ/2). 

(A9) 



The above expression can be re-written as 
l-iZ/2 l + iZ/2 l + iZ/2 l-iZ/2 



l + iZ/2 P l-iZ/2 l + iZ/2' J l-iZ/2' 

After algebraic manipulation this can be shown to be 
equal to the RHS of A4, as required. 



APPENDIX B: ELECTRO-OPTIC FEEDBACK 
VIA AN ATOM WITH JAYNES-CUMMINGS 
COUPLING AND DETUNING 

In this section we take the compound system as being 
a single mode optical cavity, with electro-optic feedback 
onto a TLA that is placed in the standing wave of the 
cavity. The Jaynes-Cummings coupling that will be used 
is V = g(au^ + era'), with g being a real constant and a 
the annihilation operator for the cavity mode. A detun- 
ing of <5<7'cr is also included. The following hierachy of 
parameters will be investigated: 



(5>.g>r>C[i? s ; 



(Bl) 



Of course, C[H S ] is really a superoperator (containing the 
system Hamiltonian terms) so here we are only referring 
to its scalar part. 

As will be shown, when the the necessary variables are 
slaved a Hamiltonian term of the form g 2 a i a/8 is ob- 
tained in the final master equation. With the above scal- 
ing, this Hamiltonian is not necessarily small compared 



to r. This makes the adiabatic elimination of the atom 
more difficult since the presumption that the atomic re- 
laxation time is much shorter than any system time scale 
is not necessarily true. To do the elimination of the atom 
rigorously we therefore transform to an interaction pic- 
ture defined by Hq = —g 2 (a i a + a'a)/8. This transfor- 
mation has the additional effect of adding a time depen- 
dence into the feedback term of the master equation. To 
nullify this we will start with a time dependent feedback 
Hamiltonian whose effect, when moved to the interaction 
picture, is time independent. The master equation in the 
Schrodinger picture is thus 

W = C[H S ]W - i[8a^a + g(aa^ + aa^),W] + TV[a]W 
+V[{aexp(-ig 2 t/8) + aUxp(igH/ 8)} a]W. (B2) 

In the interaction picture with respect to Hq the mas- 
ter equation is 



W = U*<C{H S }{UWU^)U -i g 2 {a)a + ^a)/8,W 

-i[g(aa^ + aa)) + 8^ a, W] + TV[a]W 
+V[e-xp(-iTT(T x /2)a\W. 



(B3) 

For simplification we will put A = 8 + g 2 /8. 

The expansion of Eq. (2.17) is made, with the ps now 
understood to be in the interaction picture. The time 
rates of change are 



Po 



C[H s }p + J[a]p 2 - A[a]p + Tp 2 
-ig{a?p\ 



Pl a) 



r g 2 a [ a 



■ Po 



p\ = C[H s ]px + J[a]p\ - A[a]pi + ig(p a^ 



g 2 a [ a 



-,Pi 



(B4) 

a i p 2 ) 
(B5) 



p 2 = C[H s ]p 2 + J[a]p - A[a]p2 - Tp 2 



-ig(api 



pW) 



g 2 a [ a 



■ P-2 



(B6) 



When p = po + P2 is used, the above equations give 
p = C[H s ]p + V\a\p - i [g 2 a)a/8, p] 
ig[a,px\. 



ig[a\p\] 



(B7) 



In the limit r>C [H s ] the amplitudes of p\ and p 2 re- 
spond to changes in the cavity mode much more quickly 
than po- Their equilibrium values are 



Pi 



P2 



r 

2~A 



•x(-gpa ] + g 2 {a\p 2 }), 



J[a] g 2 n t 1 



[j[a]p- ig(ap 1 - 4 fflt )} 



(B8) 



(B9) 
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These two equations can be rearranged to give pi in terms 
of pq. In the limit A ss T 2 and g 2 w TA, we find to first 
order 



P2 



Ar* 



C[a f a]) JMp+^aJMp. (BIO) 



Using this in Eq. flB8|) allows the following master equa- 
tion to be derived: 



Tq 2 

p = C[H s ]p + V[a]p + -j^V[a]p - i 



<?V) 2 « 2 



A3 



2<?* t , f, 2 9 2 r , t A * M 
C[a ( aJ ( 1 - — C[a*a\ I J[ajp. 



Ar 



(Bll) 



In the limit of A 3> T 2 , while still maintaining g 2 w 
TA, the third and fourth terms drop out, leaving the 
same master equation derived in Sec. II B, with Z — 
2g 2 a' ! a/TA of order unity. This is the same limit in 
which Walls and Milbu rn arrive at the effective Hamilto- 
nian used in Eq. (2.15) [O]. The third and fourth terms 
correspond to, respectively, an increased damping rate 
and a x^ 3 "* nonlinearity for the cavity mode. 

Note that the derived Hamiltonian term which threw 
doubt upon the adiabatic elimination process has been 
canceled. Of course, when we return to the Schrodinger 
picture it will reappea r, le aving a different master equa- 
tion from that of Sec. II B. The solution is to start with 



an extra Hamiltonian term of the form —g 2 a j{ a/ A when 
using the Jaynes-Cummings coupling. A transformation 
to the interaction picture is then not required, nor is the 
time dependence in the feedback. 



APPENDIX C: UNITARY TRANSFORMATION 

OF THE TOTAL MASTER EQUATION FOR 
ELECTRO-OPTIC FEEDBACK ONTO A MODE 

In this appendix the total master equation for electro- 
optic feedback onto a cavity is unitarily transformed so 
that the amplitude of the driven cavity may be treated 
classically, thus reducing the necessary basis size. The 
implicit form of the master equation will be used as 
this proves to be more straightforward. That is to 
say, the photocurrent will be approximated by a slightly 
smoothed version of dN/dt. Before transformation the 
implicit master equation is 



W = -i 



ga^a(b + b^ + ^{(a^ 2 -a 2 } 7 W 



-- [(-ib + ib r )I(t),W] + TV[b]W. 

We now put W = UWU^ where, 

C/ = exp[e/(t)(6-6+)/2] 



■V[a]W 
(CI) 

(C2) 



and fit) is given in Eq. ( [3 . 1 2[ ) . The unitarily transformed 
master equation is given by 



W = UU^W + WUU f + UC(U*WU)U*, 
where W = CW. Now, U is given by 

u=i(b-tf)[-rf(t)/2 + i(t)]u, 



thus the first two terms of Eq. ( C3 ) give 



-[-r/(i)/2 + /(<)] \b-b\w\. 



(C3) 



(C4) 



(C5) 



The last term of Eq. ( |C3| ) will only cause a change to 
terms that are dependent upon the driven cavity opera- 
tors. Thus, the non-linear driving and the damping of the 
system may be ignored for the present. The expression 
that needs to be simplified contains three terms (damp- 
ing, coupling, and feedback), which can be evaluated us- 
ing UbU> = ef(t)/2 + b. The damping term is 



TU 



(v[b](UiWU)\rf 



= T (r>[b]W + ef 
the coupling term is 



b-b\W 



/4 



(C6) 



-iU[g(Ja{b + tf),U^WU]U^ 
---i[gja(b + tf + ef),W], 



(C7) 



and the feedback term is 

-ib + itf)I{t),U^WU\ £7* 



iU 

T 



/(t)(6-6 t ),W 



(C8) 



Adding up the contributions from the damping, feedback, 
coupling, Eq. ( ]C5D and also the system Hamiltonian, the 
following is obtained: 



W 



ga^a{b + + ef) + ^{(a f ) 2 - a 2 }, W 



-TV[b]W, 



(C9) 



which is the master equation after transformation. There 
is now no distinction between the implicit and explicit 
forms as / is a bounded function. 



APPENDIX D: ALL-OPTICAL FEEDBACK 
ONTO AN INFINITE NUMBER OF CAVITIES 



It is of interest whether all-optical feedback can ever 
have the same effect on a system as simple electro-optic 
feedback. In this section we show that this can be 
achieved with a very large number of ancilla cavities that 
are coupled back to the system. The basic idea of the 
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all-optical feedback remains the same, in that the out- 
put of one cavity becomes the input to the next cavity. 
The cavities all have the same damping co-efficient, V. 
Damping of the system is set equal to unity. For large 
damping, the infinite number of ancilla cavities will be 
adiabatically eliminated. See Fig. [L2| 

The form of the coupling of the cavity is similar to 
that for the all-optical feedback via a single mode. It is 



Kb% 



N 



(Dl) 



where N is the total number of driven cavities, K is pro- 
portional to an Hermitian system operator and &,• is the 
annihilation operator of the j th cavity. The input and 
output fields to and from the j th cavity are represented 
by u in(j ) and u out (j). This means that u in{j+1) = u out (j)- 
The output field f rom the system is given by u ut- 

From Eq. ( |2.51 ) the contribution to the QLE for an ar- 



bitrary operator due to the shining of the — output 
field onto the j cavity can be found. The total QLE if 
there are N driven cavities is 

JV 

dr = V[J]rdt-[dVlc-c)dV ai ,r]+T^V[b]]rdt 

3=1 

JV 

-Vf^idUl^ -b]dU inU)i r] 



3=1 



JV 



3=1 



N 



(D2) 



where dll\ 



in (i) 



l in(j) 



dt. Note that the same idea as in 



Eq. ( 2.55|) has been used. The QLE for a system operator 



is 



ds = V[c*]sdt - [dV-lc - cUVin, s] 



JV 



iJ2 b ^N> s ^ dt + i ^ Hs ' s ^ dt 



3=1 



l jV 



(D3) 



Although the input fields Ui n (j) obviously depend on sys- 
tem operators, they are evaluated at a slightly earlier 
time due to the small, but finite, propagation time of the 
field from the system to the driven cavities. The system 
operator s, therefore, commutes with dUj. 

We n ow n ote that the QLE for bj has the same form 
as Eq. (HI) 



dbi 



Tb, 



VT\ 



1(3) 



iKbj 



N 



dt. 



Also, 



l in(j) 



VTbi 



(D4) 



(D5) 



To simplify matters only the non-stochastic part of bj 
and Win(j) will be considered in the derivation of the mas- 
ter equation for the system. This can be justified by 



mathematical induction. Suppose that tw^-) and bj can 
both be grouped into stochastic terms linearly depen- 
dent upo n fin and non-stochastic terms. Then it is clear 
from Eq. (D5) that u- ln /j + ij can also be group ed in such a 
manner. Therefore, in the limit in which Eq. ( D4) can be 
slaved to produce the equivalent of Eq. ( 2.54 ) it can be 



seen that bj+i will consist of non-stochastic terms, aris- 
ing from the non-stochastic terms of iti n (j) and bj, as well 
as stochastic terms linear in V{ n . To complete the math- 
ematical induction, b\ and Mj n (i) can obviously grouped 
in the manner suggested. Now, terms in bj that go as v- m 
will annihilate onto the vacuum state when the expecta- 
tion value of Eq. ( |D3| ) is taken, thus, the stochastic parts 
can be ignored as they give a zero contribution. 

An expression for non-stochastic part of Ui n (j) (denoted 
by Mj n(j) ) needs to be found in order to evaluate bj . Using 
Eq. (D5) and the slaved value bj-i it is found to be 



hn(j) 



-l + 2iK/TN\ 
l + 2iK/TN ) 



3-i 



(D6) 



Subs titu ting into Eq. ( ]D4| ) gives bj. This is then used in 
Eq. (D3). Writing Z = 4K/T, the summation term is 



l-J^il-iZ^N)- 1 



Firstly, the quotients are expanded to second order in 
Z/N. Then the contributions from the first and second 
orders are factorized, with the latter expanded to first 
order in j /N 2 . This gives 



i 

N 



j p (1 - iZ/2Ny 1 



N 



x£(l + *W (l- J -(Z/N) 2 )[Z, S ] 

3=0 



X 



(1-iZ/Nf (l- 3 -(Z/N) 



(D8) 



It is not difficult to show that the contribution of 
the j/N 2 terms is small (of order iV -1 ). Also, (1 + 
iZ/2N)~ 1 ~ 1, so we can now write the summation as 



^2j[cHl + iZ/N) j ] [Z/N, a] 

3=0 
N 

J [c f (1 + tZ/Nf] (J [(1 + iZ/N)] - I) . 



3=0 



- j 



(1 + iZ/N) 3 



3=0 



Jl^} (j [(1 + iZ/Nf 
J[J\(J[ex$(%Z)]-\)t 



Is 



(D9) 
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Terms of order TV -1 have been ignored as the limit 68,580 (1992). 

N — > oo has been taken. Returning to Eq. (D3) gives the [26] P. D. Drummond and D. F. Walls, J. Phys. A 13, 725 
same QLE as for simple feedback. In this rather imprac- (1980). 
tical way, all-optical feedback can replicate electro-optic 
feedback. 
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FIG. 1. Schematic representation of simple feedback. The 
system is taken to be a single mode optical cavity, with an- 
nihilation operator a and damping rate 7. All further figures 
will also use an optical cavity for the system. 
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FIG. 2. Schematic representation of direct detection feed- 
back onto a TLA that is coupled back to the system. The 
system damping rate has now been set equal to unity and the 
TLA damping rate is T. 
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FIG. 3. Schematic representation of direct detection feed- 
back onto an optical cavity that is coupled back to the system. 
The ancilla cavity has annihilation operator b and damping 
rate F. 
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FIG. 4. Schematic representation of all-optical feedback 
onto a TLA that is coupled back to the system. 
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FIG. 5. Schematic representation of all-optical feedback 
onto a single mode cavity that is coupled back to the sys- 
tem. A Faraday rotator (FR) and a Polarization dependent 
Beam Splitter (PBS) are included in the feedback loop. 
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FIG. 6. Wigner functions of the steady states produced 
with simple feedback for \ — n/2. Figs, (a) and (b) have 
A = 2.2, while (c) and (d) have A = 0.97. The mesh plots 
are included to aid the readers interpretation of the contour 
plots. 
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FIG. 7. Figs, (a), (b) and (c) are Wigner functions of 
the steady states produced with electro-optic feedback onto a 
TLA for x = tt/2 and A = 2.2. Fig. (a) is the adiabatically 
eliminated state. Figs, (b) and (c) represent the full dynam- 
ics with T = 2 and 20 respectively. Fig. (d) shows the Bures 
distance between the adiabatically eliminated state and the 
state produced with the full dynamics as F increases. 




X 1 Damping rate, r 

FIG. 8. Figs, (a), (b) and (c) are Wigner functions of 
the steady states produced with electro-optic feedback onto a 
mode for \ — tt/2 and A = 2.2. Fig. (a) is the adiabatically 
eliminated state. Figs, (b) and (c) represent the full dynamics 
with r = 10 and 100 respectively. Fig. (d) shows the Bures 
distance between the adiabatically eliminated state and the 
state produced with the full dynamics as T increases. The 
error bars are due to statistical error due to averaging over a 
less than infinite number of quantum trajectories. Only half 
error bars are given because, in a high-dimensional Hilbert 
space, a state with statistical errors will tend to be further 
away from the adiabatically eliminated state than the true 
ensemble average will be. 
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FIG. 9. Figs, (a), (b) and (c) are Wigner functions of the 
steady states produced with all-optical feedback onto a TLA 
for x — 7r /2 and A = 0.97. Fig. (a) is the adiabatically elim- 
inated state. Figs, (b) and (c) represent the full dynamics 
with F = 1 and 10 respectively. Fig. (d) shows the Bures 
distance between the adiabatically eliminated state and the 
state produced with the full dynamics as F increases. 
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FIG. 11. The first two plots are Wigner functions of the 
steady states produced with a x^ — t/2 non-linearity. Figs, 
(a) and (b) have driving strengths of A = 2.2 and 0.97 respec- 
tively. Fig. (c) shows the Bures distance between simple and 
"reversible" feedback for A = 0.97. Fig. (d) is a mesh plot of 
the Wigner function displayed in Fig. (a). 
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FIG. 10. Figs, (a), (b) and (c) are Wigner functions of the 
steady states produced with all-optical feedback onto a mode 
for x — and A = 0.97. Fig. (a) is the adiabatically elim- 
inated state. Figs, (b) and (c) represent the full dynamics 
with F = 1 and 10 respectively. Fig. (d) shows the Bures 
distance between the adiabatically eliminated state and the 
state produced with the full dynamics as F increases. Half 
error bars are used for the same reason as in figure H. 
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